Mutual friction in helium II: a microscopic approach 
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We develop a microscopic model of mutual friction represented by the dissipative dynamics of a 
normal fluid flow which interacts with the helical normal modes of vortices comprising a lattice 
in thermal equilibrium. Such vortices are assumed to interact with the quasiparticles forming 
the normal fluid through a pseudomomentum-conserving scattering Hamiltonian. We study the 
approach to equilibrium of the normal fluid flow for temperatures below 1 K, deriving an equation 
P_i , of motion for the quasiparticle pseudomomentum which leads to the expected form predicted by 

the HVBK equations. We obtain an expression for the mutual friction coefficient B in terms of 
microscopic parameters, which turns out to be practically independent of the vortex mass for values 
arising from diverse theories. By comparing our expression of B with previous theoretical estimates, 
we deduce interesting qualitative features about the excitation of Kelvin modes by the quasiparticle 
O ' scattering. 
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I. INTRODUCTION 



When a sufficiently fast rotating sample of liquid helium is cooled below the lambda temperature, all the rotation 
of the superfiuid becomes concentrated in a uniform array of quantized vortex filaments parallel to the axis of rotation 
[l|, |2J. By contrast, the macroscopic superfiuid velocity field, corresponding to spatial averages over regions large 
compared with the spacing between vortices, yields the usual configuration of solid body flow, v s (r) = fl 10 tz x r 
for a rotation frequency f2 ro t around the z axis. Just as the superfiuid flow is microscopically formed by vortices, 
the normal fluid consists of superfiuid quasiparticle excitations, phonons and rotons, the average flow of which is 
characterized by the normal fluid velocity field v„. In equilibrium both fluids move with the same velocity v„ = v s 
and such a behavior arises, from a microscopic viewpoint, from the vortex motion with the normal fluid velocity in 
order to avoid dissipation. That is, any relative motion of vortices with respect to the normal fluid in their vicinity, 
is subjected to a friction force that causes such a motion to eventually cease. Such a mutual friction force 3] between 
the two fluids plays then a central role in the mechanism which maintains the stability of the above equilibrium 
^__l | state. A well-known phenomenological model for this macroscopic dynamics is represented by the so-called Hall- 
■^J" • Vinen-Bekharevich-Khalatnikov (HVBK) equations [J, 0], basically consisting of a Navier-Stokes equation for the 
normal fluid and an Euler equation for the superfiuid, which, in the absence of pressure and temperature gradients, 
are coupled together only by mutual friction. There is a sinrplc configuration which allows to show the basic features 

of this process, namely rectilinear flows of uniform vorticity [5J, 
-t— > ' 

v.(r,t) = -2n.(t)y* (1.1) 

V ■ v «( r ^) = -2fi n (*)yx (1.2) 

-a- 

(y < 0), where fl s (t) and fi n (t) should converge for t — » oo to a steady state value. Then, to make contact with the 
standard rotational configuration, we may identify such a value with the former angular velocity f2 ro t . Note that this 
assignment leads to a uniform vorticity V x v s = 2 f2 rot z which coincides with that of the rotational scheme. The 
;> HVBK equations for such flows are very simple and read 



Pn -g f =F = - Ps — , (1.3) 

where p n and p s denote the normal fluid and superfiuid mass densities, respectively, and the mutual friction force F 
can be written for temperatures below 1 K as [5j, 

F = -Bp n fi rot (v„ - v s ), (1.4) 
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being B a dimensionless dissipative coefficient. Such a low temperature regime corresponds to p n <C p s , which, 
according to (|1.3[) . implies that the main time dependence should lie within the normal fluid velocity. This suggests 
that a suitable approach to the problem may consist in regarding the superfluid component as a thermal equilibrium 
heat bath which interacts with a nonequilibrium normal fluid flow. Keeping such a picture as our basic premise, we 
shall analyze in the present paper a microscopic model of mutual friction, which reproduces the main features of the 
above macroscopic dynamics, yielding an explicit expression of 5 as a function of microscopic parameters. 

The microscopic basis of mutual friction remains as one of the most intrincate problems of superfluidity. In such a 
context, theoretical approaches may strongly differ, even in significant questions such as the existence of a nondissi- 
pative component of the mutual friction force [6| , which is absent from our modelling [5] . A better understanding of 
the microscopic principles governing mutual friction would also contribute to clarify important issues on the subject 
of quantum turbulence at finite temperature [7|. In fact, it is just the mutual friction force which accounts for the 
strong locking between superfluid and normal fluid along the turbulent cascade, where recent simulations have shown 
that the residual slip velocity v s — v„ plays a central role [8J . In addition, such simulations suggest that the cross-over 
between zero-temperature and finite temperature quantum turbulence occurs at a lower temperature than the usual 
estimation of 1 K, hence partially placing the latter regime within the temperature range of the present investigation. 

Another important source of controversy arises from the mass of quantized vortices. On the one hand, many works 
have considered it as a negligible parameter under the assumption that it should be equivalent to the hydrodynamic 
mass of a core of atomic dimensions |l|. Another theories, however, yield several orders of magnitude higher values 
for the vortex mass, casting doubt on models based on massless vortices Q. Moreover, it has been argued that an 
unambiguous vortex mass may not exist, and that inertial effects in vortex dynamics may be scenario-dependent [lOj. 
Finally, we should also mention that there have been conflicting results for the vortex mass in superconductors as well 
11]. A possible way out to such uncertainties has been recently suggested based on the concept of pseudomomentum. 
In fact, just as the momentum corresponds to the generator of ordinary translations, the so-called pseudomomentum 
generates translations which keep an eventual background medium fixed [12J . Thus, the pseudomomentum often 
appears as a useful tool in fluid mechanics [13]. In the present case we shall concentrate our attention on the 
pseudomomentum of quasiparticles forming tbe normal fluid and on the pseudomomentum of vortices, being the 
background medium the superfluid. The quasiparticle pseudomomentum per unit volume turns out to be simply 
Pn(y n — v s ), which may be easily interpreted as the ordinary momentum of a normal fluid moving relatively to the 
superfluid. On the other hand, the vortex pseudomomentum requires a more careful treatment, recently shown to lead 
to an alternative approach to the dissipative vortex dynamics, free of the ambiguities carried by the uncertainty of the 
vortex mass [5| . Our present study of the dissipative normal fluid dynamics will utilize the concept of pseudomomentum 
as an important tool, finding again results which are practically independent of the vortex mass for a wide range of 
values. 

The dynamics of the normal fluid flow is far less known than that of the superfluid, because, apart from very recent 
efforts [14| , there is a lack of experimental observation techniques at these low temperatures. Only recently, theoretical 
research has been focused on this issue by means of numerical simulations, which provided valuable information about 
normal fluid flow patterns due to the mutual friction interaction |15| |. Another interesting aspect, which has been 
recently investigated, concerns the stability characteristics of the normal fluid flow under mutual friction forcing from 
the superflow [161 ]. 

Our approach will consist in assuming a heat bath formed by a vortex lattice in thermal equilibrium, which interacts 
with a quasiparticle flow. The dominant contribution to the heat capacity of such a lattice should arise from the 
thermal excitation of helical waves, corresponding to effectively independent vortices (l7| . The role of such oscillations 
in mutual friction has been scarcely treated in the literature. We are only able to mention a couple of papers [la.[l9j. 
that long ago reached the conclusion that the damping of vortex oscillations due to phonon scattering, should not 
modify appreciably the value of the friction coefficient calculated for a rigid vortex. The same conclusion was recently 
obtained for a high-frequency branch of helical waves, within a wider temperature range, including a roton-dominated 
regime [201 ]. 

In building a theory with massive vortices, one can readily make use of a close analogy with the well-known 
electrodynamical problem of a point charge subjected to magnetic and electric fields [211 ] . Particularly, the quantization 
of the theory, which greatly simplifies the treatment when the scattering excitation of vortex waves is taken into account 
[18J ] , arises immediately from this analogy. Such an analogy also leads to an immediate identification of the vortex 
pseudomomentum, allowing us to build a proper form for a pseudomomentum-conserving scattering Hamiltonian. 

This paper is organized as follows. In the next section, we propose a Hamiltonian model for the interaction of 
a vortex line with a quasiparticle gas and a background superflow of uniform vorticity. In Sec. IIIII we obtain the 
equation of motion for the quasiparticle pseudomomentum, which leads to the expected form predicted by the HVBK 
equations. In Sec. II VI we study our expression for the mutual friction coefficient B comparing this result with previous 
theoretical estimates. Finally in Sec. [V] we give our concluding remarks. 



II. HAMILTONIAN MODEL 
A. Vortex Hamiltonian 

Let us consider an otherwise rectilinear vortex filament performing helical oscillations about its unperturbed position 
parallel to the z axis [l| . The wavelength A is supposed to be much greater than the amplitude (radius of the helix) , 
and to have a full description of the helix, one should also know the direction (right or left) of the helical deformation, 
or equivalently, the direction of the wave vector fcz (k — ±2tt/A). Thus, it suffices to know the location of the vortex 
core at a given plane normal to the z axis, say z = 0, to have a full determination of the position of the whole vortex 
filament. Periodic boundary conditions over a length L (vortex line length) along the z axis determine the possible 
values of the wave vector as k = 27rm/L, where m is an integer. The vortex core position r may then be written 
as a summation over generalized two-dimensional coordinates r k associated to normal modes labeled by the wave 
vector fcz. We shall first consider true oscillatory modes (k ^ 0), since modes with k — corresponding to rigid 
displacements of the vortex filament require a separate treatment [5[ . The equation of motion for r k (t) is given by 

r k = fiz x r fe - flu k r k , (2.1) 

with 

Q = p s n/m v (2-2) 

and 

w fc = ^-[-ln(|*|o) + 0.116], (2.3) 

47T 

where m v denotes the vortex mass per unit length and k = h/m 4 denotes the quantum of circulation, given by the 
ratio of Planck's constant and the mass of one 4 He atom. The vortex core parameter a ~ 1 A in (|2.3|) is assumed to 
be much less than the wavelength (|A:|a <§C 1), which in turn ensures that to k <C fl [lil20|- The first term in equation 
(|2.ip corresponds to the Magnus force, whereas the second term stems from the induced velocity on the vortex line 
element by the helix curvature. Assuming that the vortex has a counterclockwise circulation, the frequencies (|2.2j) 
and (|2.3|) are positive, and then equation (|2.1[) turns out to be analogous to that ruling the two-dimensional motion of 
a negative point charge, in a uniform magnetic field parallel to the z axis and subjected to a harmonic central force. 
Such an equation derives from the following Hamiltonian: 



^(v 2 fc + ^ fc r£), (2.4) 



v fe = — + -zxr fc , (2.5) 



being 



where pfe denotes the conjugate momentum to r^. In fact, from Hamilton equations it is easy to check that (|2.5[) 
corresponds to the velocity rfe, while the acceleration is indeed given by (|2.ip . From the above electromagnetic analogy 



it is also useful to represent the coordinate r k as the sum of the center coordinate Rq of the cyclotron circle plus 
the relative coordinate R fc from such a center [21| . 

r fc = R fc) + R' fc = R Q fc) + v fc x z/O. (2.6) 

(k) 

Then, from Hamilton equations one easily obtains the following pair of coupled equations for Rq and R fc : 

R fe) = -u k 2 x H { k) -uj k zx R' fc (2.7a) 

R' fe = (Q + uo k ) z x R' fe + uj k z x R { k) . (2.7b) 

The solution of the above system can more simply be expressed in complex notation (r = x + iy) , and to first order 
in uj k /Q, we have, 

R^ = Al exp(-Mdkt) - ^-A 2 cxp[in(l+io k /n)t} (2.8a) 

R' fe = A 2 exp[itt(l + uj k /n)t]- ^Aiexp(— iu k t), (2.8b) 



where A\ and A2 are complex numbers depending on initial conditions. Then, neglecting corrections of order uj k /H,, 
the solution corresponds to a decoupling of the system (|2.7[) : 



R« = - Wfc zxR fc) 



Rl 



n z x r'. 



A" 

,(fe) 



(2.9a) 
(2.9b) 



that is, we shall simply have circular trajectories, namely for R a clockwise one with angular frequency uj k and for 
R' fc a counterclockwise one with angular frequency fl. Then, it becomes clear that Rq should be ascribed to Kelvin 
modes, while R' fc should correspond to the cyclotron ones [E1201- Theoretical estimates of the vortex mass [l|, Q lead 
to cyclotron frequency values in (|2.2[) of the order or greater than ksT/h for T < 1 K. This seems to indicate that 
quantum effects could be of importance. Quantization of coordinate and momentum arises straightforwardly from 
the electromagnetic analogy and reads 21(, 



R 



(fe) _ 



Rj, 



Pa 



2p s nL 



2p s nL 
lhp s K 



-, — ik 



ifc *a fe (x + *y)+e^4(x-*y)] 



-«*(/£- a fc )(*- 



*yj 



Ahz 



(4-/3 fe )(x- Jy )] 



(2.10) 

(2.11) 
(2.12) 



where a, (0l) denotes a creation operator of right (left) circular quanta. The z dependence of coordinates and 
momentum corresponds to the rotation generated on following the helix path. Normal modes corresponding to fcz are 
ruled by a Hamiltonian of the form (I2.4[) : 



H, 



Q u k r 



(2.13) 



dz ^r[v k 

where, taking into account (|2.10|) to (|2.12|) . one may verify that vf and r| do not depend on z and also that the limit 
io k /n<^ 1 yields g(|, 

'/,"'■ ~ ' '' "~i '"/."i t- 77)- (2-14) 



H k = hn{a\a k + -) + hw k {plli K y 



That is, both polarizations (cyclotron and Kelvin modes) become decoupled, as seen from a classical viewpoint in 
equations (J2.9I) . Here it is instructive to evaluate the limit of a massless vortex, which is often found in the literature 
(l7J |. In fact, the limit Q —)■ 00 in (j2.ll) transforms such an equation into a first-order one, with the consequence 
that the Cartesian components x k and y k of r^ become conjugate variables. Note that this amounts to ignoring 



the cyclotron motion by setting r^ = Rq ' (cf. Eq. (J2.6I) ). where the components of R in (|2.10[) obey canonical 
commutation relations. 

The Hamiltonian of the k — modes reads Q 



Ho 



m„L 



>0a 



u 0y) 



£l rot p s K,Ly 



with vq x = pox/m-v — Qyo, vq v = po y /m v . The quantization of the coordinate rp = R ( 



(0) 



Rn — R 



(0) 



(2.15) 
vo x z/fl 



remains given through the previous expressions f|2. 10[) and (J2.11I) , whereas the momentum expression (|2.12[) becomes 
changed according to the Landau gauge as p = yJhp s n/(2L) [i(/3g — /3o)x + (aj + ao)y]. The Hamiltonian (|2.15|) 
can be exactly solved [21] yielding in the limit J7 rot <C Q the decoupling of cyclotron and translational modes: 

"■"rot , 



H = Ml(ala + -) 



: (^-/3o) 



(2.16) 



2' 2 

Thus, the Heisenberg equation of motion for the cyclotron coordinate R is given again by (|2.9bp . whereas the 
translational coordinate evolves according to the superfluid velocity field 



j(0) 



R 



Z \ £ ro ^ 1 q X 



In conclusion, the vortex Hamiltonian may be written as follows: 

r 



h v =j: 



Ml 



a{a k 



f^k (fitPk + \ 



Ml r 



-{Pl-M 



(2.17) 



(2.18) 



with k = 2ixmjL and m an integer. 



B. Quasiparticle Hamiltonian 

The normal fluid will be represented by the following Hamiltonian: 

# n = ^/!w q a q a q , (2.19) 

q 

where a], denotes a creation operator of quasiparticle excitations of pseudomomentum ftq and frequency w q . Such a 
frequency corresponds to the lab frame and may be written as a Doppler-shifted frequency from the superfluid frame, 
w q = Uq + q- v s , where to q is the familiar (isotropic) dispersion relationship of 4 He quasiparticle excitations. Now, 
the background superflow velocity should be much less than the Landau critical velocity ~60 m/s, so we may safely 
approximate w q = uj q in (|2.19p . Note also that we disregard any interaction between the quasiparticles themselves, 
since we shall work at low enough temperature, so that they remain dilute allowing their treatment as a noninteracting 
gas. 

C. Interaction Hamiltonian 

The interaction Hamiltonian between the vortex and the quasiparticles will be represented by the pseudomomentum- 
conserving form: 

H^t =5Z / dz A pq fl pOqexp[-i(p - q) ■ r(z) - i(p z - q z )z], (2.20) 

P,q J ° 

whe re r( z) = R (z) + R'(z), being R (z) = £ fc R fc) (2) and R/(z) = £fc R fe( z ); with R { Q k) (z) and R' k (z) given 
by (|2.10p and (|2.11|) . respectively. The parameters A pq in (|2.20p represent scattering amplitudes depending on 
wave vectors of scattered quasiparticles. The vortex pseudomomentum per unit length [5] K(z) = —p s nz x Ro(z) 
integrated along the vortex line yields the generator of vortex translations or vortex pseudomomentum L K.(z)dz = 

—p s KLz x Rg , which involves only translational coordinates, as expected. Then, adding such a pseudomomentum 
to the quasiparticle pseudomomentum ^ ^q a q aq, we have the pseudomomentum of the whole system, which can 

be shown to commute with H lnt . Note that only the x-component of the vortex pseudomomentum p s kLYq will 
commute with the vortex Hamiltonian (|2.18p . unless J7 rot = 0. This result may be easily interpreted, since a superflow 
of velocity v s — — 2 f7 rot y x produces a translation symmetry breaking in the y-direction. 

III. EQUATION OF MOTION FOR THE NORMAL FLUID FLOW 

The interaction Hamiltonian (J2.20I) is difficult to deal with, so recalling the low amplitude of the helical oscillations, 
we may rewrite the exponential in (|2.20p as 

cxpH(P - q) • r ( z ) - i(Pz - Qz)z] = 

exp[-i(p - q) • r[, 0) ] exp[-i(p z - q z )z] exp[-i(p - q) ■ (r(z) - r[, 0) )] (3.1) 

and next approximate the last exponential to first order in r(z) — R . This procedure, however, is not valid for vortex 
modes with frequencies approaching zero, i.e. the lowest part of Kelvin's spectrum, as noted early by Fetter |l8[. In 
fact, he showed that retaining a finite number of terms of such an exponential expansion leads to divergent results, 
analogous to those of the "infrared catastrophe" in electrodynamics. Here it is expedient to recall that within our study, 
each vortex forms part of a vortex lattice which will be regarded as a heat bath for the quasiparticle flow. Now, it is well 
known that rather simple models of heat bath often provide suitable descriptions of realistic environments [22| . Relying 
on this hypothesis and to overcome the above difficulty, we shall represent Kelvin's spectrum by a single frequency Wq, 
which will be eventually regarded as a temperature-dependent parameter in order to take into account the distinct 
features of the interaction of such waves with phonons and rotons. In summary, we shall make use of a simplified 
model of heat bath consisting of vortex modes of two frequencies (wq <C Q) of opposite polarization. On the other 
hand, neglecting the vortex displacement in the y-direction ,5|, we shall replace the translational coordinate operator 

in the first exponential on the right-hand side of Eq. (|3.ip by the c-number (R ') = — 2 O rot (Yq )tx (cf. Eq. (|2.17p ). 
Such a replacement becomes equivalent to having time-dependent scattering amplitudes in Eq. (|2.20p . i.e. with a 



time dependent phase factor, A pq exp[i(p x — q x )2fl lot (F )i\. This factor, however, would not have any practical 
incidence, since our results will be shown to be dependent upon the absolute value of the scattering amplitudes. 
Finally, taking into account these approximations we may replace exp[— i(p — q) • r(z)] ~ 1 — i(p — q) • (r(z) — Rq ') 
in Eg.' (|2~20f yielding (20} 



Hint = J -£— ^2 8 P *qApcia i p a cl {[(qy-py) + i(q x -p x )}[al + {l-5ko)l3k] 
V ^,p,q 

+ [(p»-g»)+*(?x-p*)][ak + (l-<7fco)#i]}. (3-2) 

There is an additional parameter to be taken into account in our vortex heat bath, that is the total number of 
modes 2L/A m i n , where the factor 2 in front of the expression comes from the two possible signs for k. We shall assume 
for simplicity that both polarizations have a common 'ultraviolet' cutoff A m i n , which according to Sec. Ill Al should 
be greater than the vortex core parameter (~ 1 A) and the mean radius of the helix. Such a radius turns out to be 
of the order of the core parameter for cyclotron modes, while for Kelvin modes in a lattice of f2 rot ~ 1 s" 1 has been 
estimated [17| as ~ 10 3 A y/T/K. So we shall assume A m ; n ~ 10 3 A in our calculations. 

A total Hamiltonian given by the sum of (|2.18|) . (|2.19[) and (|3.2[) yields dissipative evolutions for each vortex mode, 
which for cyclotron modes are ruled by [20| : 

(Hi) =n(i-—)zx <R4) - n—(K), (3.3) 

V PskJ p s K 

where explicit expressions for the longitudinal and transverse friction coefficients, D and D', respectively, can be 
found in [23| . There we have shown that a scattering amplitude given by [24[ 



Ap^-^^KIKI, (3-4) 

where V denotes the volume of the system and ui' (c s ) denotes the guasiparticle group (sound) velocity, leads to a 
very good agreement with the experimental determinations of the longitudinal friction coefficient for temperatures 
below 1 K. The presence of the viscous force represented by the last term in (|3.3|l leads to a vanishing asymptotic 
velocity, (R' k } t ^,oo = 0, and the same behavior will present the coordinate, (R,j.)t_ >0O = 0. 

To summarize, we have focused in previous works on the dynamics of an individual vortex line which interacts with 
a guasiparticle heat bath. Now, being focused on the time evolution of the normal fluid, we shall concentrate on 
the dissipative dynamics of a guasiparticle flow which interacts with the heat bath formed by a uniform array of N v 
guantized vortex filaments. We have derived in Appendix [A] the system (|A6[) of non-Markovian eguations that rule 
the time evolution of the guasiparticle populations n q . Then, from (JA6I) one easily obtains the following eguation of 
motion for the guasiparticle pseudomomentum: 

V hqhq = —?- V |A pq | 2 <5 p ^(p-q)(p-q) 2 / dr cos[(uj p - u q + Wj)r]{n q (i - r) 

q PsKAmin p q . J 

x [H-n p (t-r)] + [n q (t-r)-n p (t-r)][e' ta "/ fcsT -l]- 1 }, (3.5) 

where Wi denotes the frequencies f2 and wq. We shall assume that the guasiparticle numbers in the above expression 
are well described by a local equilibrium form: 



,/fe B T _ ,1-1 , fiq-[vn-v a ] 



q 4k B Tsmh 2 (hu; q /2kBT) 

where the quasiparticle freguency is measured from a reference frame where the local normal fluid is at rest, w q = uj q — 
q- [v n — v s ]. Thus, the guasiparticle pseudomomentum J? fiqn q = ALp n (y n ~ v s ) corresponds to a macroscopically 
small area A of the x-y plane containing N v vortices, where the spatial dependence of the fields v n and v s can 
be neglected. Then, assuming a time dependence stemming exclusively from v n (i) (cf. Sec. HJ, a straightforward 
calculation leads to the following non-Markovian equation: 

t 
d7-[v„(t-T)-v.]A*(T), (3-7) 



II 



with a memory kernel given by 

k B Tp n p s n 2 X 



mO) = , T r °2 X Yl l A pql 2(S P,?, cos [( w P- w 9 + u '0' r ]{|p-q| 4 5+K) 



mm 

p, q, « 

|4 ia .. /„ .. ~\i4 



where 



+ [|z x (p x z)| 4 - |z x (q x z)| 4 ] g-(wi)} n(u q ) [1 + n(w„)] [1 + ra^)], (3.8) 

n(w p ) , n(w g ) exp[7i(w g - w p - Wi)/k B T] 

9±{wi) = —, — r± ; , x (3.9) 

n(uj q - Wi) n{io p + Wi) 

and n(w) — [exp(fiw / k B T) — 1] — . In the thermodynamic limit, the summations over p and q in Q3.8P become 
integrals and /z(t) acquires a hnite lifetime. If such a lifetime can be regarded as microscopic in comparison with 
the observational timescale, equation (I3.7[) may be transformed according to the Markov approximation into the 
differential equation: 

v„ = -v\\ n {t) - \,\ (3.10) 

with 

V = / / u(r)dT = — — ^- V |A p q| 2 6 p ,qJ(u p - Ug + Wi) 

JO K B 1 PnPsK Amin 



p, q, 



x |p-q| 4 nK)[l+n(w p )][l + nK)], (3-11) 

where the continuum limit corresponds to the replacement X) P q $pzQz ~~ ^ [A 2 L/(2ir) 5 ] Jd 3 p/d 3 q S(p z — q z ). Actu- 
ally, we have studied in Appendix [B] the non-Markovian equation (I3.7[) finding that memory effects are negligible for 
f^rot <C Wo, which will be assumed hereafter. Finally, taking into account (|3.10p . (|1.3[) and (|1.4[) . we may obtain the 
expression of the mutual friction parameter from B = i//Q mt . 

IV. STUDY OF THE MUTUAL FRICTION PARAMETER B 

An explicit expression for the friction parameter B can be extracted by computing the right-hand side of equation 
(j3~TTj) (see Appendix©: 

19ft 3 



B = 70^) 2 c 2 p n p s k B TX mm ? ^ + "Ml I * Kl^ + «*> t 1 + »M E T <*> ^ (4 ' X) 



where 



L[P > q) -\q 2 p(p 4 + q 4 /5 + 2p 2 q 2 ) (q < p) [ ^ Z) 

and q = q) denote the roots of the equation ui q = u> p + wi. From (j4.ll) we may see that B consists of two terms 
arising from the frequencies Wi — wq and Wi = O. Such contributions, however, are weighted by respective factors 
[1 + n(wo)] 2> [1 + u(O)], so the cyclotron contribution will be always negligible with respect to that arising from 
the frequency wq and we have that, in practice, B will correspond to the limit of massless vortices, fl — > oo. The 
expression (|4.1j) leads to simpler phonon and roton approximations. If we restrict ourselves to wq < 10 10 s _1 , such a 
frequency may be neglected everywhere in (j4.ll) . except in the factor [1 + n(wo)]. Then, to approximate for phonon 
temperatures (T<0.4 K), we use the linear dispersion relation u p = c s p and get 

254.9 [l + n(w )](k B T) 3 

ph = fi 2 r 4 n \ . ' ( 4 - 3 ) 

n c s p s A mm 

On the other hand, for temperatures above 0.6 K, only the portion of the dispersion curve around the roton minimum 
makes a significant contribution to the integrand in (|4.1j) . then making use of the usual approximations in roton 
calculations [2J], we obtain 

2.079 [l + n(w )]hk$y/l^T 
1 = c 2 /Tip X ■ ' ( } 
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where n and fco are parameters entering the Landau parabolic approximation, uj p — A/H+ H(p — ko) 2 /2/z. 

Experimental determinations of B have been reported only above 1.3 K. However, we may utilize the following 
expression valid for temperatures below 1 K [l| , 

B = ™ (4.5) 

and replace D in (J4.5I) by means of the Iordanskii theoretical estimate for the phonon temperature range [25| , yielding 

B ph = 8.17 -^. (4.6) 

7714 Cg 

On the other hand, for the roton temperature range, we may replace D — p n VG<J\\ in (|4.5[) yielding, 



ft^JM-iJiVfK-l, (4.7) 

K V 7T/J 

being vq the average group velocity of rotons and cni ~ 8.38 A the roton scattering length [l|, @|. Notice that we have 
replaced the friction coefficient D by expressions corresponding to a straight vortex, since corrections due to vortex 
bending should be negligible, as seen in Sec. HI This of course does not mean that vortices would remain straight 
against the quasiparticle scattering; on the contrary, thermal excitation of Kelvin waves is undoubtedly expected to 
occur, although details of this process are not evident from expressions (|4.6[) and (|4.7[) . However, some features about 
such a process may be deduced from our results (|4.3[) and (J4.4I) . First it is convenient to discuss the physical meaning of 
the frequency w . Since such a frequency is intended for representing the whole spectrum of Kelvin waves (|2.3|) in the 
context of an interaction with quasiparticles, it should not be surprising to find that quite distinct values of wq could 
be required in order to better estimate interactions with phonons or rotons. This amounts to assuming a dependence 
of wq on temperature, which may be fully extracted by equating our results with the theoretical expressions (|4.6[) and 
(|4.7|) . In fact, from (14.31) and (|4.6|l we may conclude that phonon scattering at a temperature T should be expected 
to excite Kelvin waves about a representative frequency given by 

w ^ 5.64 x 10 8 s~ 1 K~ 3 T 3 , (4.8) 

while (|4.4|) and (J4.7I) imply that roton scattering should be likely to excite Kelvin waves about frequency 

w„ ~ 2.08 x lO^s^K -1 T. (4.9) 



Recall that according to the Markov approximation (Sec. IIIip one should assume wq 3> f2 ro t ~ 1 s , which sets 
up a lower bound for the validity of the result (|4.8|) at temperatures above ~ 0.01 K. In addition, the assumption 
Amin ~ 10 3 A implies an upper bound for the Kelvin spectrum, max(ojfc) ~ 10 s, which turns out to be consistent 
with a phonon temperature range below 0.4 K in (|4.8p . However, the values arising from (|4.9p seem to be overestimated 
for roton temperatures, since they would only be consistent with a A m i n of order 10 2 A. In addition, such values of 
Wo could reach the order of cyclotron frequencies arising from some theories of the vortex mass 9] , contradicting the 
assumption uj^ -C Q of Sec. Ill Al This suggests that only the qualitative trend wq ~ T should be taken into account 
from the result (|4.9[) . 

V. CONCLUDING REMARKS 

We have analyzed a microscopic model of mutual friction represented by the dissipative dynamics of a normal fluid 
flow, which interacts with the helical normal modes of vortices comprising a lattice in thermal equilibrium. Such 
vortices interact with the quasiparticles forming the normal fluid through a pseudomomentum-conserving scattering 
Hamiltonian. Assuming a simplified model for the vortex heat bath, we have derived an equation of motion for the 
quasiparticle pseudomomentum leading to the expected form predicted by the HVBK equations. We have shown 
that the mutual friction coefficient B turns out to be practically independent of the values of vortex mass arising 
from diverse theories. Finally, from a comparison of our expression of B with previous theoretical estimates, we have 
deduced interesting qualitative features about the interaction of quasiparticles with Kelvin modes, namely phonon 
(roton) scattering at a temperature T should be expected to excite Kelvin waves about representative frequencies 
proportional to T 3 (T). 
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Appendix A: Derivation of the equation of motion for the quasiparticle populations 

Our starting point is the Heisenberg equation for the quasiparticle number operator a^ a q : 



^K a q) = v/ry— H S P*qA A P<l[(Px-<lx)+i(<ly-Py)} a l a ci [<4 3)t + (! - M/3i j) ] 
' s k,p,j 

+ A pq [iPx ~ &) - i{ly ~ Py)\ l a k ] + (1 - 4o)/3£ j)t ] aq a P + A Pq [(P* ~ 3x) ~ *(«* ~ Pa)l a p a q [ a !fc 

+ (1 - 4o)/3i i)f ] + A pq [(p a - q x ) + i(q y - Pj/ )] [<4 i)f + (1 - *«,)/#>] ot a p }, (Al) 

where j labels each vortex (1 < j < N v ) and we have taken into account that only the interaction Hamiltonian (|3.2[) 
yields a nonvanishing commutator with such an operator. Next we write out the following expressions for the operators 
appearing on the right-hand side of the above equation, which arise from the formal solutions of the corresponding 
Heisenberg equations: 



a p (i)a q (i) - e^-^) t a p (0) aq (0) + l W-^- f dre^~^ £ 6 P ^ A q , p {[fo, - q' y ) +i(p x - q' x )] 

V PsK J ° k.q'.j 

„U)U+ — ^ _i_ n _ X, \aU)f + _ ^M I \(J _ „ \ , ;r m _ J \i \JS)(+ _ _\ 1 /1 _ £. _ , \fl(i)t * 



x [<#"(* - r) + (1 - 4o)tf J (t - r)] + [((?; - jfc) + ife - <£)] [ajf j (* - r) + (1 - 6 M )fi>\t - r)]} 

— 111 —1— 1 1 1 I — #1.1 ■ ' ( ' " ' 

+ [(Qv 1y) + i(q' x &)] [<*?>(* - r) + (1 - M/f )f (* - r)]>4 (t - r) a q ,(i - r), (A2) 



x oL(t - r) a q (t - r) - <^ A qq , {[(</ - g„) + i(<£ - q x )] [a[ j) \t - r) + (1 - 4o)/?i J) (i - r)] 



„0')t 



(t) = e iat aj^O) + JJt- J dre«» £ S^A^ [{p' x q' x ) + t(p£ - q' y )\ a p ,(< - r^C* - r) (A3) 
V 5 ° P ', q ' 



(A4) 



(3^{t) = e^'fiWQ) + \U— k J dre^ £ 5 p , q , A p , q , [& - £) - i(p' y - <^)]o^(t - r)<v(i - r 

V S ° p',q' 

The above expressions and their Hermitian conjugates are then replaced on the right-hand side of (|A1[) . retaining only 
second-order terms in the scattering amplitudes (weak-coupling approximation). The resulting expression becomes 
greatly simplified when taking its ensemble average according to the following prescriptions: 

(i) All vortex-quasiparticle correlations are neglected, i.e. any average of a product of vortex and quasiparticle 
operators is approximated by the corresponding product of vortex and quasiparticle separate averages. Such a 
procedure may be regarded as arising from, (a) an assumption of vanishing initial correlations, i.e. assuming a whole 
system density operator given by a product of a vortex operator and a quasiparticle operator, and (b) the above 
weak-coupling approximation by which such vortex-quasiparticle averages should be calculated to zeroth-order in the 
scattering amplitudes. 

(ii) The assumption of a vortex heat bath corresponds to vortex operator averages represented by thermal equilib- 
rium expressions, i.e. (4 i)ta fc j) ) = [exp(m/k B T)-l]- 1 , (/3[ m i3^ ] ) = [exp^wo/fcsT)-!]- 1 and (a[ m ) = (^f )f ) =0. 

(iii) We assume that the system is close to equilibrium so that all nondiagonal quasiparticle averages can be 
neglected, i.e. 

(a P <z q ) = S pci n p . (A5) 

Thus, we arrive at the following system of non-Markovian equations for the quasiparticle populations n q : 

-?T = i r-H |A Pq | 2 <W(p-q) 2 {/ drcos[(«,-« p + «; i )T]{np(t-r)[l + n q (t-r)] 

at n p s KA m in Jo 

P ,t 

+ Mt-^-n^-T)]^'/** 7 "-!]- 1 } 

drcos[(w p -u: q + Wi )t] {n q (i - r)[l + n p (t - r)] + [n q (t - r) - n p (t - r )][ e W« T - lp 1 }}, (A6) 
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where it;, denotes the frequencies ft and wq. It is easy to verify that (|A6|) fulfills ^ driq/dt — 0, in agreement with 
the quasiparticle number-conserving Hamiltonian. It is also instructive to rewrite (|A6j) in the Markovian limit, 



dn^ = 2nL 2 N v ^ |A pq | 2 * p . & (p - q) a {5(«, -«, + t«0 {np(t)[l + «,(«)] + [np(t) -«,(«)] [e to «/ kBr - l]" 1 } 

at np s KA mm — 

- tf(wp - uj q + m,) {ra q (i)[l + np(t)] + [n q (i) - n p {t)][e hw *' kBT - l]" 1 }}, (A7) 

and verify that thermal equilibrium populations, n p = \e hu) p/ k BT _ i]- 1 ^ i eac { to a vanishing result on the right-hand 
side of the above equation. 

Appendix B: Study of memory effects 

To explore the validity of the Markov approximation, it is convenient to Laplace-transform equation (|3.7p according 
to v(z) = J e lzt [v n (t) — v s ]dt, (Im z > 0). Then one easily finds that 

v(z) = [v n (0)-v a }/\fi(z)-iz), (Bl) 

where the Laplace transform of the memory kernel can be written as a Cauchy integral: 

2m J_^ u-z 



with 

2A min A:_BTp„/9 s K 



r/(w) = r °* ^ |A pq | 2 ^g i8 5(a; p -u;g + uJiH-u;){|p-q| 4 5 + (uJi) + [|zx (pxz)f 



p, q, * 

AM4 



- |z x (q x z)| 4 ] 3_(wi)} n(w,)[l + n(w p )] [1 + n(w,)]. (B3) 

The poles of v(z) in (|B1[) arise from the equation z = —ifi(z) and the Markov approximation corresponds to the 
solution Zm = —iji(z — > i0 + ) = —iv [cf. p. lip ]. A non-Markovian solution can be found iteratively, i.e. we may 
begin with the Markovian ansatz zq = —ip,(0) and then proceed with z\ — —ifl(zo), Z2 = —iji(zi), and so on. Then 
the solution z s is better worked out in terms of the Taylor expansion of p,(z) around the origin: 

z s = -iu[l - »j5'(0) - /i"(0)A(0)/2 +•••], (B4) 

where the second and third terms inside the square brackets represent first- and second-order corrections to the 
Markov approximation, respectively. The Cauchy integral (IB2|) and its derivatives in equation (|B4|) can be written as 

M 

/t(0) = 1/(0) (B5a) 

A'(0) = ^/ ™[u(u) + u(,-w)-2v(p)] (B5b) 

2m J_ 00 w l 

jS"(0) - i/'(0), (B5c) 

where the calculation of v(oj) from expression (|B3|) may be reduced to a single one-dimensional integral (see [241]): 

19S 3 n mt 



u(u>) = 



x 

3 



140(2 7 r) 2c ^ nPs fc s TA- £ [1 + " W] J M < 1 " K + «* + ^ [1 + n( " p)] 

X) r(p,gf )g+(w i ) + ^(p,qf ) )</_(«*)> (B6) 



where <7± (wi ) is given by equation (|3.9[) with 

w q = o; p + ^i + w, (B7) 

^'W-i o 2 p(p 4 + a 4 /5 + 2pV) («<J>) ( j 
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$r^-^M-« 4 +p73 + fpV) (p<«) mq x 

*(P,?)-<{ (? 2 p(p 4_ (? 4 / 3_| p 2 (Z 2 ) (g < p) (B9) 



and g = qf denote the roots of the equation (|B7I) . The above expressions allow a direct computation of the equations 



but it will be more instructive to perform the following dimensional analysis whose results were numerically 
verified. Firstly, from equations (IB5a[) and p. Ill) we have [1(0) = v = f2 ro t-B which, according to the experimental 
values of B, must be of order J7 rot at most. Then from equation (|B5b|) it is not difficult to realize that p.'(0) should be 
~ ^rot/wo & n d similarly /i"(0) ~ Q TO t/wQ. Thus, the first- and second-order corrections in equation (JB4I) turn out to 
be, respectively, of order J7 rot /wo and (r2 rot /u>o) 2 which tells us that the Markovian limit corresponds to J7 rot <C Wq. 
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